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1. 

I.    Introduction  and  Summary . 

In  a  collection  of  evaporating  drops,  such  as  an  aero- 
sol, each  drop  is  affected  by  the  presence  of  the  other  drops. 
The  evaporation  of  these  other  drops  tends  to  raise  the  vapor 
pressure  or  concentration  in  the  neighborhood  of  any  drop  above 
what  it  would  be  if  that  one  drop  evaporated  alone.   In  effect, 
the  drops  compete  for  space  into  which  to  evaporate.   If  there 
are  N  drops  per  unit  volume  of  aerosol  then  each  drop  has  an 
effective  volume  of  N   into  which  to  evaporate.   Therefore  we 
take,  as  a  model  of  a  collection  of  evaporating  drops,  a  single 
drop  within  an  impermeable  container  of  volume  N~  .   For  mathe- 
matical convenience  we  assume  that  the  container  is  a  sphere  with 
the  drop  -  also  a  sphere  -  at  its  center.   Then  we  have  to  deter- 
mine the  radius  R(t)  of  the  drop  as  a  function  of  time.   To  do 
so  we  must  simultaneously  consider  the  concentration  c(r,  t)  in 
the  space  between  the  drop  and  the  container,  the  radius  of  which 
we  call  M. 

The  physical  principles  on  which  this  determination  is 
based  are  : 

a)  The  flow  of  vapor  is  assumed  to  satisfy  Pick's  lqw  and  there- 
fore the  concentration  satisfies  the  diffusion  equation. 

b)  The  concentration  is  assumed  to  have  a  constant  value  -  the 
saturation  value  -  at  the  drop  surface. 

c)  The  rate  at  which  the  drop  radius  decreases  is  assumed  to  be 
proportional  to  the  rate  at  which  vapor  diffuses  away  from  the 
drop  surface,  (i.e.  conservation  of  matter). 

d)  No  vapor  flows  through  the  container. 
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e)  The  concentration  is  assumed  to  be  initially  uniform  through- 
out the  space  between  the  drop  and  the  container. 

f)  Temperature  changes  and  heat  flow  are  ignored  -  the  process 
is  assumed  to  be  isothermal. 

Mathematically  these  conditions,  from  which  R(t)  and  c(r,t) 
can  be  determined,  are,  as  given  in  [6] 

(1)  c   +  -  c   =  c.     for  M  >  r  >  max[0,  R(t)],   t  >  0 
'  rr   r  r      t 

(2)  c[R(t),  t]  =  1     for     0  <  t  <  tQ 

(3)  acr[R(t),  t]   =  R(t)   for     0  <  t  <  tQ 
(Ij.)  c(r,  0)     =0     for     1  <  r  <  M 

(5)  c„(M,  t)     =0     for     t  >  0 

(6)  R(0)   =  1 

(7)  cr(0,  t)     =0     for     t  >  tQ 

In  equations  (1)  -  (7)  all  quantities  are  in  the  dimen- 
sionless  form  given  in  [3]  (see  eqs.  (7)  -  (12)  of  [3]).   There- 
fore since  the  initial  drop  radius  is  the  unit  of  length  and  the 
container  radius  must  exceed  the  drop  radius,  we  have  the  restric- 
tion M  >  1.   The  constant  t   is  the  time  required  for  complete 
evaporation,  if  it  occurs,  and  is  infinite  if  complete  evaporation 
does  nat  occur.   The  constant  a  is  a  given  positive  number  deter- 
mined by  the  coefficient  of  diffusion,  the  density  of  the  liquid 
in  the  drop,  the  initial  concentration,  the  saturation  ooncentra- 
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tion  and  the  initial  drop  radius. 

The  problem  formulated  above  reduces,  when  M  =co,  to  the 
case  of  a  drop  in  infinite  space  with  no  container,  which  was 
treated  in  [3]  and  [5],   In  [3]  approximate  results  were  obtained 
for  R(t)  and  for  t  by  formal  calculation.   In  [5]  it  was  proved 
that  the  problem  with  M  =oo  has  a  unique  solution  and  that  the 
drop  completely  evaporates  in  a  finite  time.   A  process  for  con- 
structing the  solution  was  presented  and  upper  and  lower  bounds 
for  R(t)  and  for  t  were  obtained  which  agree  closely  with  the 
approximate  results  of  [3].   For  practical  applications  the  upper 
and  lower  bounds  are  so  close  together  that  they  determine  R(t) 
and  t   with  satisfactory  accuracy. 

In  the  present  case  in  which  the  drop  is  enclosed  within 
a  container  of  finite  radius  M  the  problem  formulated  above  also 
has  a  unique  solution  for  R(t)  and  c(r,  t).   The  radius  R(t)  is 
a  monotonically  decreasing  function  of  t  but  the  drop  does  not 
alwavs  evaporate  completely.  Whether  or  not  this  happens  depends 
upon  the  container  radius  M.   There  is  a  critical  value  of  the 
container  radius  M  =  (a"1+  1)  '•*    such  that  for  M  >  M  the  drop  does 
evaporate  completely  in  a  finite  time  and  the  concentration  every- 
where within  the  container  approaches  the  constant  value 

c    =  (M/M)^   as  t  — >  oo.   When  M  =  M  the  drop  radius  R(t)  ap~ 
oo      ' 

proaches  zero  as  t  — >oo,  and  the  concentration  approaches  unity. 
If  M  <  M  the  drop  does  not  evaporate  completely,  but  instead  its 
radius  approaches  the  value  [1  -  v{W -   1)]  '^   and  the  concentration 
approaches  unity  as  t  — •>  oo.   All  of  these  results, which  are  quite 
reasonable  and  almost  obvious  physically,  have  been  proved  and  are 
reported  on  in  this  report.   However  the  existence  of  the  solution 
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is  proved  only  when  a  or  H  satisfy  certain  conditions.   Undoubtedly 
the  solution  exists  for  all  a  >  0  and  all  Ir  >  1,  but  other  methods 
are  necessary  to  prove  it. 

In  ad'ition  to  the  proofs  of  the  above  statements,  the 
main  results  of  the  present  report  pertain  to  the  determination 
of  the  radius-time  curve  R(t)  and  the  evaporation  time  t  •  It 
is  shown  that  the  function  R(t)  is  the  limit  of  the  sequence  of 
functions  P^-it)  constructed  in  section  IV.  More  useful  than  this 
fact,  however,  is  the  result  that  the  functions  p  (t)  also  provide 

upper  and  lower  bounds  for  R(t)  and  for  the  rate  of  evaporation 

• 

R(t).   Explicitly  this  result  is  that  for  any  value  of  n 

(8)  p2n(t)   <  R(t)   <  P2n+1(t) 

(9)  p2n(t)  <  R(t)   <  P2n+1(t) 

By  means  of  these  equations  arbitrarily  close  bounds  on  R  and  R 
can  be  obtained.   These  results  have  only  been  proved  for  restricted 
values  of  a  or  of  M  -  for  other  values  of  these  constants  they  have 
been  proved  for  only  a  limited  period  of  time.   As  an  example  of 
the  use  of  these  results  the  following  lower  bound  on  R(t)  has  been 
obtained  using  only  the  first  term  in  the  sequence  p  (t).   The 
bound,  which  holds  if  M  >  3,  is 

(10)  R2(t,  >1+te[A  +  ^ilog(l-^li^)] 

A  similarly  obtained  lower  bound  on  R(t)  is  given  in  equation  (68). 

Other  bounds  on  R(t)  and  on  tQ,  when  tQ  is  finite,  are 
given  in  the  following  theorem,  which  is  proved  by  means  of  Green's 
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theorem  in  section  VI. 

(11)        Theorem  8. 

a)  R2(t)      <      1  -  f~  (1  -  c      )    =      <Kt) 

v    '      —  1+a  oo '   —     Y      ' 

b)  R2(t)      >     X  -  a(M2-l)    -  2af  =    XU) 

c)  R2(t)      <     IMii^g       =        *(,}      for  4(t)    >  0 

d)  R3(t)      >     1   -   a(M3-l)    =  P3^  ,        for  M  <  K 

e)  tQ  >     ^  -  |(M2-  1)  ,        for  M  >  M 

OO 

The  expression  R  occurring  in  the  lower  bound  in  (d)  is 
the  exact  final  radius  for  the  case  M  <  M,  while  the  upper  bounds 
in  (a)  and  (c)  reduce,  in  this  case  to  the  trivial  bound  R  <  1. 
For  M  >  M,  R(t)  vanishes  at  a  finite  time  t   and  then  all  bounds 
except  (d)  yield  significant  information. 

In  Figures  1,  2,  and  3  th9  applicable  bounds  on  R(t)  are 
shown  for  the  cases  M  =20,  60  and  200  respectively,  i.e. 
container  radius  =20,  60  or  200  initial  drop  radii.   The  constant 
a  has  been  chosen  as  a  =  2  x  10   ,  corresponding  to  the  evaporation 
of  a  drop  of  water  in  air.   For  this  choice  of  a,  M  =  36.8I4.,  so 
that  in  the  case  M  =  20  <  M  the  drop  does  not  evaporate  completely. 

In  this  case  (Fig.  1.)  the  drop  radius  is  bounded  from 
below  by  the  constant  R   .   We  also  have  the  lower  bounds  ty(t) 
(good  initially)  and  X(t)  (  a  transition  bound).   The  upper  bound 
4=1  merely  states  that  the  drop  evaporates  while  P(t)  is  a  fair 
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bound  during  its  short  period  of  applicability. 

In  Pig.  2,  for  M  =  60  >  M  ,  the  drop  evaporates  in  finite 
time  (2.32  x  lcr  <  t  <  3.22  x  10  *")  and  again  Theorem  7  supplies 
rough  bounds  over  the  full  range,  while  the  iteration  gives  \|/(t) 
an  excellent  lower  bound  initially  but  applicable  for  only  a  short 
time.   The  curve  K  (t),  (included  for  comparison),  represents  the 
actual  solution  for  M  =  oo  since  the  upuer  and  lower  bounds  for  this 
case,  (see  [5],  P«  2)  coincide  within  the  accuracy  of  this  graph. 

For  M  =  200,  Pig.  3,  the  upper  bound  <K  t )  coincides  (with- 
in graphical  accuracy  with  (X(t)  while  now  it  is  the  lower  bound 
\|((t)  which  yields  a  better  result  than  X  (M  • 

It  seems  likely  that  much  better  bounds  can  be  obtained 
than  those  given  here,  and  therefore  the  most  important  results 
of  this  report  are  the  methods  and  not  the  actual  bounds.   In 
fact,  the  same  methods  have  been  applied  to  the  corresponding  one 
dimensional  problem( evaporation  of  a  slab  in  a  bounded  region)  and 
the  bounds  obtained  were  excellent.   It  is  to  be  hoped  that  simi- 
lar results  will  be  obtained  in  the  present  problem  as  well. 

II.    Reformulation  of  the  Problem. 

Instead  of  c(r,  t)  it  is  more  convenient  to  deal  with  the 
function  u(r,  t)  defined  by 

(12)  u(r,  t)    =   rc(r,  t)   . 

Prom  (l)-(-7)  we  obtain  the  following  equivalent  set  of 
equations  for  u(r,  t)  and  R(t) 
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(13)  urr  =   ufc    ,        for  n  >  r  >  max[0,  R(t)  ,   t  >  0 

(11+)  u(R(t),t)  =  R(t)  ,     for  0  <  t  <  t  . 

(15)  aur(R(t),t)  =  R(t)R(t)  +  a  ,   for   0  <  t  <  tQ 

(16)  u(r,  0)   =  0     ,    for  1  <  r  <  M 

(17)  Mur(M,  t)  -  u(M,  t)  =  0       for   0  <  t 

(18)  R(0)   =   l 

(19)  lira  —   exists  and  is  bounded  for  t  >  t 
r — >0  r 

ru  -u 

(20)  lim  — =■£-  =   0  ,      t  >  t 

r 

In  the  range  0  <  t  <  t   the  conditions  (19)  and  (20)  can  be  dis- 
regarded.  We  now  seek  two  functions  u(r,  t)  and  R(t)  which  sa- 
tisfy (13)-(20). 

III.   Method  of  Solution 

The  problem  formulated  in  (13-lS)  is  a  free  boundary  value 
problem  for  the  heat  equation.   To  solve  it  we  follow  the  general 
method  for  solving  such  problems  given  in  [ly] . 

The  first  step  in  this  method  is  to  construct  an  auxili- 
ary function  up(r,  t)  which  we  now  proceed  to  define.   Let  p(t) 
be  a  smooth  function  such  that  p(0)  =  1.   We  use  the  notation 


up(p-,  t)  =  lira   up(p(t)  +  e,  t) 

0<e — >0 

uP(p±,  t)   =    lira   u£(p(t)  +  e,  t) 
r  0<e—>0 
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The  function  up(r,  t)  is  then  defined  for   -  co  <  r  <  M,  r^p(t), 
t  >  0  by  the  following  conditions: 

(21)  u£r     =     up  r^p(t),      -oo<r<M 

(22)  up(p+,    t)    -  uP(p-,    t)      =      p(t) 

(23)  a[up(p+,    t)    -  uj(p",    t)]    =     p(t)p(t)   +  a 

(2k)  up(r,   0)      =     0    ,      r  <  M    ;        lira         up(r,    t)     =     0 

r — >  -  cd 

(25)  Mup(M,    t)    -   up(M,    t)      =      0    ,         t    >   0      . 

To  find  up,  we  write  it  as  a  sum  vp  +  wp,  where  vp  is 
defined  for  -  co  <  r  <  oo  and  satisfies  equations  (21)  to  (2l\.) , 
while  wp  is  a  compensating  part  defined  for  r  <  M  in  such  a  way 
that  the  sum  also  satisfies  the  condition  (25).   The  unique 
solution  for  vp  has  been  found  in  [5],  page  7,  equation  (19).  As 
to  wp,  we  try  to  represent  it  in  the  form 

t 

(26)  wp(r,  t)   =   -  f     [M  vp(M,  tf)  -  v(M,  tf)]a(r,  t  -  <T)dcT  . 

i 

Here  a(r,  t  -  d")  is  to  be  determined  in  such  a  way  that  wp  satis- 
fies (21)  and  (2i|)  while  vp  +  wp  is  to  satisfy  (25).   We  then  get 
the  requirements  a   =  a . ,  and  a(r,  0)  =  0  to  satisfy  the  former 
condition,  while  to  satisfy  the  latter,  Mar(r,  t  -  6)   -  a(r,  t  -  o*) 
must  be  a  "delta  function".   Since  (Mar  -  a)  also  satisfies  the 
heat  equation,  we  now  find,  by  using  the  well  known  properties  of 
the  fundamental  solution  of  the  heat  equation,  see  [1],  that 
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(27)  Mar(r,    t   -  6)   -  a(r,    t   -  6)   =  - 


1        r-M 


'-M   \ 


?  (t-<5)^  ^2./^ 


where 


(28) 


ty(x)      =-•      n 


1/2   exp    (-x2)    . 


Solving    (27)    for   a,    we   get,    on  using  the   condition   a   — >  0      as 
IH    — >    oo    , 


a(r,    t  -  6)  =  -»y  exp( 


S*^/ 


r-M v/t  -cf 

2/tTo 


^(d)do' 


-  00 


Inserting   this   in   (26)    and   adding   it   to   vp  we   now   get,    after 
some   simplifications, 


(29)    up(r,    t)   =  r 


z(r,t) 


oo 


f  tj(tf)dtf  +    /  Y,(d)dd 

z(r,0)  ?(r,0) 


-A 


h(z(r,0))-nU(r,0))J 


r 

.  |_  r  £(<aj?iia  ^j  +  ^und* 

J  i/t-T 


0 


-   2(M-l)exp(^  -  2/151^)     ^  n«f)d< 

+   £(r,0) 


ou 
M 


2_ 

aM 


r 


00 


/"    p«J)p(rf)exp(i=|  -  ^  *)(    /  n(T)dT)dcf 

0  M  ^ 


. 


' 
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Here 


r  z  =    z(r,C)    =   r"P(tf) 
2/E^cT 


(30) 


z(r,0)    = 


V 


2/t-cf 

+    oo      if      r  >   p(t) 

z(r,t)   =  o        if     r  =   p(t) 

-   oo      if     r  <   p(t) 


r-1 
2/t 

2M-r~l 

2/t 


It  is  observed  that  the  formula  (29)  for  up  is  defined  for  r  =  p(t), 
and  we  shall  say  that  up(r,t)  takes  on  this  value  for  r  =  p(t). 
Differentiating  (29)  we  get 


(31)  uP(r,t)  =_y 


z(r,t) 


h(c)dcr  5-  /  h(o-)dd  - 

z(r,0)  ?(r,0) 

t 


1 
2/t 


rj(z(r,0))-  ty£(r,0)] 


o 

00 

2(M-1)  a    /  t    2/t^(r,0)\  p       ^i*\a* 


M 


+  S(r,0) 


_2 
aM' 


t  oo 

f    P(c5)p(0)exp(£=2  -  ^4-^  *)(  X— ¥T)dT)d(f  ' 


Formula  (31)  for  up  has,  likewise,  a  meaning  when  r  =  p(t),  and  we 
shall  say  up  takes  on  this  value  when  r  =  p(t)  . 

Concerning  the  function  up  we  now  state  the  following 
Theorem  1 

a)  up,  up  are  continuous  for  r  <  p(t),  and  p(t)  <  r  <  M,  t  >  0 

b)  up  is  the  unique  solution  of  problem  (21)  through  (2$)    . 
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c)    While  up  and  u£  are  discontinuous  for  r  =  p(t),  we  have, 

u?(p±,  t)   =  uP(p,  t)   +  |(t) 

u£(p±,  t)   =  UP(p,  t)   +  ^(a^pttJpCt)  +  1)  . 

This  theorem,  as  well  as  ather  theorems  not  proved  here,  are  proved 
in  [10]  . 

The  remaining  theorems  assert  that  u  with  p  replaced  by  an 
appropriate  R(t)  represents  the  solution  to  the  original  problem. 
R(t)  is  the  evaporation  curve.   It  is  defined  as  the  solution  of 
the  equation 


(32)       Mu£(p",  t)  -  uP(p",  t)   =  0  ,   p(0)  =  1 


Theorem  2 

Equation  (32)  has  a  unique  solution  R(t)  at  least  in  the  inter- 
val 0  <  t  <  min  (t  ,  T)  where  t  is  such  that  R(t  )  =0  and  T  is 
defined  in  Section  IV  by  equation  (55)*      In  this  range  R(t)  is  con- 
tinuous, while  it  is  differentlable  in  the  open  interval,  and  is 

•         -1/2 
monotonically  decreasing;  near  t  =  0,  R(t)  =  0(t  '    )  . 


Theorem  3 » 

uR(r,    t)   =  0        for  r  <  R(t),        0  <  t   <  min(tQ,   T)    . 
Theorem  \\% 

u(r,    t)    =  uR(r,    t)        for  R(t)    <r<M,      0<t<  min(to,T) 


12. 
satisfies  equations  (13)  through  (18)  and  hence,  with  R(t)  deter- 
mined by  equation  (32)  constitutes  the  unique  solution  of  the  prob- 
lem in  the  indicated  rantTe , 

It  will  be  shewn  later  that  both  possibilities  t   >  T  and 

c  o 

t   <  T  can  arise,  depending  on  the  value  of  M.   Thus  for  that 
range  of  M  for  which  t  <  T  we  have  obtained  the  complete  solu- 
tion, since  we  can  extend  u  for  t  >  t   by  the  same  process  as  that 
used  in  [5]«   When  t   >  T  we  have  obtained  an  incomplete  solution 
of  the  mathematical  problem  although  the  results  may  still  be 
sufficient  for  any  practical  applications. 

In  proving  Theorems  1  -  k   we  find  that  Theorem  1  is  almost  a 
direct  consequence  of  a  corresponding  theorem  in  [5]    although  cer- 
tain modifications  are  necessary.   It  is  in  attempting  to  prove 
Theorem  2  by  tne  iteration  method  that  we  are  forced  to  restrict 
the  time  interval  (see  Sec.  iv-2).   Theorem  3  requires  an  entirely 
different  proof  from  that  used  for  the  corresponding  theorem  in 
[5]»   It  is  necessary  to  apply  the  Strong  Maximum  Principle  of 
Nirenberg  [8]  and  also  a  theorem  suggested  by  Nirenberg  concerning 
the  normal  derivatives  of  solutions  of  parabolic  equations  near  a 
maximum  point  which  occurs  on  the  boundary.   The  theorem  is  analo- 
gous to  that  of  Fopf  [2]  for  elliptic  equations. 

The  proof  of  Theorem  !|.  follows  the  lines  of  the  correspond- 
ing proof  in  [5]  • 


13. 

IV.    The  Equation  for  R(t)  and  its  Solution 

The  equation  for  the  evaporation  curve, that  is  equation  (32), 
Involves  certain  limits  which  can  be  evaluated  by  using  Theorem  lc. 
We  then  get  instead  of  (32), 

(32<)  M(u£(p,t)  -  ^   pp  -  \)    -  (up(p,  t)  -  |p)  =  0,   p(0)  =  1 

This  may  be  rewritten  in  the  following  form 


(33)     pp  =  V[p]  =  2a|  u£(p,  t)  -  i-  uP(p,  t) 


+ 


hn> 


p(0)  =  1 


The  functional  V[p],  is  explicitly,  on  using  (29)  and  (31), 


(31*) 


-co 


oo 


drf 


V[p]    =  2a(l   -   4^J)[    /      ^|(tf)d^+     y     ty<0 

z(p,0)  S(p,0) 


-   a(-i-  -  ^)     [h(z(p,    0)    -  *l(S(p,0)) 


/  ^^  [  (-^w-(«^H 


do-. 


At  this  point  it  is  possible  to  explain  why  equation  (32)  was 
chosen  for  the  determination  of  R(t)  while  according  to  the  general 
theory  of  [l\.]    other  equivalent  equations, apparently  simpler,  could 
have  been  used,  and  while  in  the  case  of  the  unconfined  drop  dis- 
cussed in  [5]  the  equation  u£(p~,  t)  =  0  was  used  successfully 
(see  equation  (22)  there).   If  one  would  follow  the  routine  of  [5], 
the  equation  for  R(t)  would  be  much  more  complicated,  since  all  the 


Ill- 
exponential  factors  appearing  in  equation  (31)  would  be  reflected 
in  it.   It  is  only  by  using  combinations  of  terms  appearing  in 
(32)  that  one  manages  to  get  rid  of  a  number  of  unwieldy  exponen- 
tial terms. 

In  solving  (33)  we  are  guided,  to  some  extent,  by  the  tech- 
niques used  in  [5].   In  [5],  Kolodner  produces  an  alternating 
sequence  of  functions,  which  converges  to  the  solution  P(t),  by 
exploiting  the  fact  that  the  operator  appearing  in  his  equation 
is  a  "decreasing  functional'1.   His  method  has  been  extended  (see 
[7])  to  include  operators  which  can  be  written  as  the  sum  of  a 
"decreasing11  and  an  "increasing"'  functional.   In  our  case,  for  the 
operator  V[p]  in  (3^4-)  we  can  write 


(35) 


V[p]      =     D[p]    +   I[p] 


where 


(36)      D[p]    = 


-2a    [      AtycOdtf-    y\(<5)ddl    -  £-^(z(p,0))-  i^(?(p,0) 


/ 


/ 

0 


p(<5)p((5)     \    zh(z) 
/t-o*  (   /t-cr 


1 

M" 


/n  J     v/t-o     /z-6    J 


and 


15. 


(37)  i[p]  =  ^MLiL 


00  00 

A}(tf)dtf  -  f  Y|(d)d(5 
-2(p,0)    S(p,0) 


+  -£-tyS(p,0))  +  ^^(z(p,0)) 
/t 


0 


'p(c<)p(tf)  \  _  _X__ 
i/t  -  cf   (   /E"^ 


i*jg>-¥«l- *&-*«]}«■ 


1'he  class  C  of  functions  p  for  which  D  and  I   are  "decreasing'1  and 
"increasing''' ,  respectively,  is  constructed  in  [10],  where  the  follow- 
ing theorem  is  proved. 
Theorem  5» 

If  p-,,  p2  e  C  for  0  <  t  <  t   and  if 


*   (1) 
Pi  Pi  1      P2  P2 


then 


Pi  <  P2 


Pi  <  P2 


and 


D[P]_]   >  D[p2]   ,     ItPil   1  I[p23 


for  0  <  t  <  t  . 

The  class  C  is  composed  (essentially)  of  all  continuously  dif- 
ferentiate positive  decreasing  functions  p(t)  for  which 


(1)   Throughout  this  paper,  when  comparing  two  functions  f(t)  and  g(t), 
the  statement  f  <  g  shall  mean  f(t)  <  g(t)  for  each  t  in  the  specified 


interval. 


16, 


(38)  z(P,  tf)  >  -  1//?   . 


The  restriction  (38)  is  necessary  even  in  [5]  and  is  practically 
equivalent  to  requiring 


(39)  p(t)p(t)  >  -  p(l  +  -M 

/nt 


with 


(14-0)  P  •   - .38 

"    i+yB7i 

In  attempting  to  set  up  an  iteration  scheme  the  key  problem  is 
that  of  determining  a  lower  bound,  that  is  a  function  W(t)  e  C  such 
that 

(la)  D[l]  +  I[W]   >  WW 

Then  'or  any  p1 ,  p2  e  C,  if  p1p1  <  0  and  p2p2  >  WW  we  would  have, 
by  Theorem  5 

(1+2)    D[p]_]  +  I[p2]   >  D[l]  +  I[W]  >  VJW   . 

We  were  not  able  to  construct  a  W(t)  such  that  (l+l)  holds  for  all  t, 
and  this  is  the  cause  of  the  time  restriction  discussed  above.   The 
V.'(t)  that  we  shall  use  is  constructed  as  follows. 
On  assuming  0  <  W(t)  <  1,  we  can  show  that 


V 


-         -       .       -    ■  I   ■  ' 


■■}■     i    . 


17. 


(k3)        D[l]  +  I[W]  >  -  A(t)  +  /  W(cf)W(tf)  .  f(t,  6) 


dd 


c 


where  the  following  definitions  have  been  used. 

M-l 


ihh) 


A 


(t)   =  a(z  /     h(d)dd  +  — — )      , 
^  0  ^ 


(kS) 


(U6) 


0   ^ 


—  f(t,  d)dd  =  -i- 


^  •* 


1  -  B(t) 


B,t)  -  i-a^f     ¥*w*-^(#j-™.) 


0 


If 


(1+7) 


i/tf  W(tf)W((5)   >  /tW(t)W(t) 


then   (I+3)    implies 


(lj.8)      D[l]    +   I[W]   >   -  A(t)    +  N/t  W(t)W(t)      /     --  f(t,d)dd 

l/5 


0 


or 


(lj.9)      D[l]    +   I[W]      >     -  A(t)    +  W(t)W(t)    (1  -  B(t))    . 


Therefore,  if  we  define  '■  ( t )  by 


' 


. 


18, 


W(t)W(t)   = 


A(t) 
"  BTtJ 


(5o) 


w(o)  =  1  , 


with  A(t)  and  B(t)  given  by  (i\k)    and  (i|6)  we  see  immediately  upon 
substitution  in  (1+9)  that  (1+2)  is  satisfied.   It  can  also  be  veri- 
fied that  (l|7)  is  valid  for  this  choice  of  V  . 

p 

Unfortunately,  the  function  B(t)  has  a  zero,  t_,(t  '^M  /3 ) 

and  hence  WW  is  defined  only  for  0  <  t  <  t   . 

We  now  define  the  sequence  of  functions  f  p  (t)rand  the  sequences 
5a  )  andfk  >  of  numbers  as  follows  : 


(5D 


P0(t)p0(t)  =  W(t)W(t) 


0  <  t  <  a. 


Po(0)   =   1 


(52) 


px(t)  =    i 


For  n  =  1,  2,  . .  • 


DtP2n-l]  +  I[P2n-2^   '   °  *   t   *   a2n-l 


/2n+l(t)p2n+l^) 


.  a2n-l  -  t  -  a2n 


P2n(0)   =   1 


D[p2n]    +   I[P2n-l]  >      0   <   t    <   a2n+1 


(Sk)       X^l(t,p2n+l(t>a,    c 


>      a2n+l  1  t  1  T 


P2n+1(°) 


=   1 


19. 

Here  a,  =  min(aT,  T),  where  al"  is  the  smallest  zero  of  W(t)  and 
T(<  t  )  Is  defined  as  the  unique  positive  solution  of 

(55)  Wtl)W(T)   =  w(T)w(T)   ; 
w(t)  in  (SS)    is  defined  by 

(56)  w(t)w(t)   =   -  p(l  +  -J-  )   , 

(57)  P  -  /gL    . 

This  definition  of  p  is  dictated  by  the  method  of  proof  of  one  of 
the  theorems.   Also,  (57)  places  a  further  restriction  on  a.   We 
must  have  (see  ( i+0 )  )  p  <  (1  +/57ii)~  ,  i.e. 


(58)  a     — - P  ~  .095 


For  n  =  1,  2,  ... 


a2n  =  min(a*n,  T)   ,   a2n+1  =  min(a*n+1  ,  T) 


where  a~    is  the  smallest  zero  of  p2n(t)  and  a2n+l  is  the  smallest 
zero  of  P2n+1(t)p2n+1(t)   . 
Also 

(59)  ko  =  °    *     kl  =  1 


' 


20. 
(60)         k2n  =   P2n(a2n-l)    ■    k  2n+l  =   p2n+l(a2n+l> 

Most  of  the  above  definitions  are  illustrated  in  the  schematic 
drawing  (Pig.  U.)  . 

We  now  have  the  following 
Theorem  6. 

For  a,  p  satisfying  (57),  (58)  there  exists  a  positive  number 
T  defined  by  (55)  such  that 

a)  The  sequence   a   converges  either  to  T  or  to  a  number  t   <  T 

b)  The  sequence  k   converges  either  to  zero  (when  t   <  T)  or 
to  a  number  K  >  0  where  R(T)  =  K  . 

c)  For  any  t  <  min(t  ,  T),  there  exists  an  N  such  that  for 

n  >  N  each  p   is  defined  in  0  <  t  <  t  and  the  sequence  p  converges 
to  a  function  R(t)  which  is  continuously  different! able  in  the  range 
0  <  t  <  min(t  ,  T).   R(t)  is  the  unique  solution  of  (8l)  and  R(t)  <  0  . 

d)  For  any  n 

p2n  -  p2n+2  -  R  -   ?2n+3   -  p2n+l 

•       •         •      •         • 

P2n  -   p2n+2  -  R  -  p2n+3  -  P2n+1 

where  the  inequality  between  any  two  functions  holds  throughout  their 
common  domain  of  definition. 


.  ■ 


' 


I 


21. 

V.    Estimates  for  R(t)  From  the  Iteration  Method. 

Theorem  6  supplies  us  with  a  sequence  of  bounds  for  R  and  R.< 
However  these  bounds  are  extremely  complicated  in  form  and  hence  ws 
shall  now  produce  some  more  usable  bounds  by  simplifying  a  few  of  the 

p  .   We  have 
"n 

(61)       R(t)R(t)   >  p0(t)po(t)   =  W(t)W(t)   =  -  ||*|   , 
It  is  easily  shown  that 


(62)       A(t)  <  a(l  +  -i-  .)   , 

/nt 

while 


(63)     B(t)  >  l  -^  -^yt  > 


These  results  hold  if 


(6k)  M   >  3 


and 


««>        *  <  i(^)2  <  *B  • 


On  integrating  (61)  we  obtain 


l(W2(t)  _!)>./   c(l+l//*cr)    dJ 

0    1-  — ^|/o 

U       M-l    M70 


22 


or 


(66)  R2(t)     -  W2U)     ■  1   ♦  !±SM  E(t)      =     vj;2(  t ) 

,/i 

In    (61+)    E(t)    is   defined   by 


(67)        E(t)=  *  +  Sfiy10^1-^^^  • 

Prom  this  we  obtain 

(68)         R(t)   >   -|f|j(l  +  ^E(t))-1/2   . 


V 


71 


These  bounds  were  referred  to  in  the  Introduction. 

a 

In  order  to   obtain  upper   bounds    for  R  and  R  we    shall   estimate 
p->    •      By   definition 

C  D[p2]    +   I[p1]  ,        0   <   t   <   a3 

33  /  0  ,        a     f  t   <   T 

but,    in   the   range    0  <   t   <  a-,    ,    since    ppp2   >  ww   and   P^   =   1>   we   have 

* 

P3P3     1     Dtw^    +   If1^      ■ 

On  eliminating  some  of  the  small  negative  terms  and  making  some 
further  simplifications  we  obtain  the  following  upper  bound: 


1 


• 


23. 


(69 


[erfc    (1/2/F)    -  erfc    (^tl) 


1   -   erfc 


("#)] 


Vt 


} 


where 


(70) 
(71) 


f(t)     = 


^ + " > 


w(t)     =     (l  -  U/t  f(t))1/2 


and 


(72) 


2        *  ?        U) 

erfc(x)    =   1   -   erf(x)    =   1   -  —      T         exp(-d    )dcf 


In  order  to  obtain  upper  bounds  on  R  and  R  from  (69)  it  is  necessary 

•  2    2**  7~2~ 

to  integrate  PP  numerically,  to  get  R  <  P   and  then  R  <  PP//P  . 


(1) 


These  are  tabulated  functions.   See  e.g.  Pierce  [9]. 


■ 


2fc. 

VI.    Estimates  for  R(t)  from  Green's  Theorem. 

By  again  applying  Nirenberg's  theorems  (see  Section  II)  we 
can  establish  the  following  theorem  concerning  the  solution  c(r,  t) 
Theorem  7» 

a)  0      <     c(r,    t)      <      1 

b)  ct(M,    t)      >     0    . 

1  ,    for   1   <  M  <  M 


c)      lim  c(r,    t)    =  c 


oo 


fc->°°  Km/m^/3 


(M/M)  /J,  for  M  >  M 


where 


(73)  m  =  (^)1/3   • 

(For  a  =  2  x  10~  ,   M  =  36.814).   Physically,  part  a)  states  that 
the  vapor  density  remains  between  the  initial  density  and  the  sa- 
turation density.   Part  b)  states  that  the  vapor  density  at  the 
container  wall  increases  with  time  as  the  drop  evaporates.   Part 
c)  gives  the  exact  final  value  of  the  density  for  both  incomplete 
and  complete  evaporation.   Theorem  7  holds  wherever  the  solution 
c(r,  t)  exists . for  all  t  >  0  . 

We  shall  now  establish  a  number  of  bounds  for  R(t)  and  for  the 
evaporation  time,  t  ,  (where  it  exists).   The  bounds  are  contained 
in  Theorem  8  which  is  stated  in  the  Introduction. 
Proof  of  Theorem  8. 

We  apply  a  special  case  of  Green's  Theorem,  (see  e.g.  [1]): 
for  any  regular  solution  u(r,  t)  of  the  equation  u   =  u.  in  a 


. 


- 


. 


.  t     I 


■    - 


:  '  '       . 


. 


i 


- 


closed  domain  D  we  have,  for  any  constants  a  and  b, 

(7k)  0  I   [ar  +  b]u  dr  +  [ ( ar  +  b)  ur  -  aujdtj  =   0  . 

On  substituting  the  data  of  our  problem  into  formula  (7k)    we  obtain, 
(after  simplifying  and  stating  in  terms  of  c(r,  t) 


(75)  SiMf[l  -  R3(t)]  ♦  |[i  -  R2(t)]| 

t  M 

=  bt  -  b  J  c(M,t)dT  +   /   (ar  +  b).r.c(r,  t)dr 
0  R(t) 

All  parts  of  Theorem  8  now  follow  from  (75)  by  choosing  appropriate 
values  for  a  and  b  and  applying  Theorem  7.  For  example,  let  a  =  0, 
b  =  1  in  (75)  to  get 

t  M 

(76)  5+i[l  -  R2(t)]  =  t  -  J    c(M,T)dT  +   J  rc(r,  t)dr  . 

0  R(t) 

Since  c(M,  T )  <  c    ,   c  >  0  ,  we  can  write 

(77)  ^i[l  -  R2(t)]   >  t  -  C(X)t   . 


From  (77)  part  a)  follows.   Parts  b),  c),  d)  are  proved  by  similar 
arguments.   Parts  e)  and  f)  are  proved  by  setting  R( t  )  =0  in  (75) 
and  then  applying  the  estimate  s. 

By  using  bound  f )  it  can  be  shown  that  for  M  >  275,  T  >  t   so 
that  in  this  case  the  iteration  procedure  of  Section  IV  gives  the 
complete  solution. 


.    ■ 


26. 


It  should  be  noted  here  that  although  the  method  of  this  sec- 
tion produces  some  good  bounds  for  R(t)  and  t   it  gives  no  infor- 
mation  whatever  on  the  rate  of  evaporation  R(t).   Thus  the  results 
of  Section  V  supply  the  only  bounds  on  R(t). 


27. 
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